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ABSTRACT Preeclampsia is one of the most severe gestational complications which is one of the leading causes 
of maternal and perinatal morbidity and mortality. A growth in the incidence of severe and combined forms 
of the pathology has been observed in recent years. According to modern concepts, inadequate cytotrophoblast 
invasion into the spiral arteries of the uterus and development of the ischemia-reperfusion syndrome in the 
placental tissue play the leading role in the development of preeclampsia, which is characterized by multiple- 
organ failure. In this regard, our work was aimed at studying the patterns of placental tissue transcriptome that 
are specific to females with PE and with physiological pregnancy, as well as identifying the potential promising 
biomarkers and molecular mechanisms of this pathology. We have identified 63 genes whose expression proved 
to differ significantly in the placental tissue of females with PE and with physiological pregnancy. A cluster of 
differentially expressed genes (DEG) whose expression level is increased in patients with preeclampsia includes 
not only the known candidate genes that have been identified in many other genome-wide studies (e.g., LEP, 
BHLHB2, SIGLEC6, RDH13, BCL6), but also new genes (ANKRD37, SYDE1, CYBA, ITGB2, etc.), which can be 
considered as new biological markers of preeclampsia and are of further interest. The results of a functional an- 
notation of DEG show that the development of preeclampsia may be related to a stress response, immune proc- 
esses, the regulation of cell-cell interactions, intracellular signaling cascades, etc. In addition, the features of the 
differential gene expression depending on preeclampsia severity were revealed. We have found evidence of the 
important role of the molecular mechanisms responsible for the failure of immunological tolerance and initia- 
tion of the pro -inflammatory cascade in the development of severe preeclampsia. The results obtained elaborate 
the concept of the pathophysiology of preeclampsia and contain the information necessary to work out measures 
for targeted therapy of this disease. 

KEYWORDS microarrays; placenta; genome-wide analysis; preeclampsia; transcriptome; gene expression. 
ABBREVIATIONS DEG - differentially expressed genes; MFDs - multifactorial diseases; PE — preeclampsia; 
GWAS - genome -wide association study; SNP - single nucleotide polymorphism. 



INTRODUCTION 

The numerous genome-wide association studies 
(GWAS) conducted so far have provided valuable in- 
formation on the genetic architecture of multifactorial 
diseases (MFDs) and revealed hundreds of the risk alle- 
les of single-nucleotide polymorphisms (SNPs) associat- 
ed with many phenotypes. However, they explain only 
a relatively small part of the inheritance of complex 



traits and have only a very mild impact on the pheno- 
type of associated variants [1]. These results raise the 
missing heritability problem, which is being intensively 
discussed today. Another limitation of the GWAS effec- 
tiveness related to studies of the hereditary component 
of predisposition to MFDs is associated with the use of 
tagSNP. The risk alleles identified in GWAS typically 
do not belong to the "causal" ones, but are in linkage 
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disequilibrium (LD) with functionally significant vari- 
ant alleles [2]; therefore, the biological interpretation of 
GWAS results is a serious problem. 

The current approaches to the identification of the 
"causal" allelic variants linked to the polymorphisms 
detected in GWAS are based on the analysis of the cod- 
ing or transcribed genomic regions [2-4]. However, the 
vast majority of SNPs identified in GWAS are located 
in the non-transcribed regions. They are not linked to 
variants located in exons, and the mechanism of their 
action is apparently associated with the regulation of 
gene expression [5, 6]. Therefore, post-genomic meth- 
ods (which readily provide information on almost all the 
components coordinating the basic functions of genes, 
RNA, and proteins at different hierarchical levels) be- 
come especially relevant in studying the genetic archi- 
tecture and molecular mechanisms of MFDs. One such 
approach, namely the high-performance measurement 
of gene expression using microarray technology, was 
used in the present work to characterize the transcrip- 
tome patterns in normal pregnancy and preeclampsia 
(PE), one of the most severe gestational complications. 

Preeclampsia, which is associated with the multiple 
organ dysfunction syndrome, is a specific syndrome 
that occurs after the 20 th week of pregnancy and is 
characterized by hypertension and proteinuria. PE is 
diagnosed in 70% of hypertensive disorders in preg- 
nancy, and an increase in the incidence rate of severe 
and combined forms of this disease has been observed 
in recent years [7]. Despite the large number of theories 
related to etiopathogenesis (neurogenic, hormonal, pla- 
cental, immunological, genetic, etc.), numerous studies 
of the mechanisms of development of this disease, and 
the emergence of new therapies, PE remains a leading 
cause of maternal and perinatal morbidity and mortal- 
ity. The disease is responsible for up to 70% of stillbirths 
and miscarriages; the risk of perinatal losses increases 
almost fivefold in PE [7, 8]. 

According to the modern concepts, the etiopathogen- 
esis of preeclampsia is closely related to inadequate cy- 
totrophoblast invasion in the uterine spiral arteries and 
development of the ischemia-reperfusion syndrome, 
which induces oxidative stress and systemic inflam- 
mation [9, 10]. Etiological factors and the mechanisms 
of this disorder remain unclear and require close atten- 
tion. In order to identify the likely candidate biomark- 
ers of PE and to study the molecular mechanisms of 
gestational complications, we analyzed the patterns 
of the placental transcriptome that are specific to PE 
and physiological pregnancy, since the placental tis- 
sue obviously plays the key role in the development of 
PE. The strategy of using microarrays in this context 
seems to be reasonable and powerful enough, as it al- 
lows one to thoroughly investigate the possible changes 



in gene expression associated with the pathophysiology 
of preeclampsia at the transcriptome level. 

EXPERIMENTAL 

Characteristics of the examined groups 

A total of 10 patients with PE and 11 patients with 
physiological pregnancy (the control group) were ex- 
amined {Table 1). The questionnaire included demo- 
graphic information (ethnicity) and anthropometric 
parameters (height, weight), lifestyle (smoking habit, 
psychoactive substance abuse), as well as information 
about the somatic and obstetric-gynecological history. 
PE was diagnosed based on leading clinical symptoms 
of various severity, such as proteinuria, edema, hyper- 
tension (systolic blood pressure above 140 mm Hg, di- 
astolic blood pressure above 90 mm Hg) according to 
the 10 th revision of the International Classification of 
Diseases. PE severity was evaluated according to the 
criteria of the clinical protocol 2012 "Hypertension dur- 
ing Pregnancy. Preeclampsia. Eclampsia" [11]. 

The group of PE patients was heterogeneous both 
in terms of severity (the study included six patients 
with moderate and four patients with severe PE) and 
the presence of prior diseases and comorbidities. Four 
patients were diagnosed with PE in the absence of 
background diseases; in others the gestational compli- 
cations developed against the backdrop of extrageni- 
tal diseases including hypo-/hypertensive type neu- 
rocirculatory dystonia, chronic pyelonephritis, chronic 
cholecystitis, and chronic arterial hypertension. Six 
females in the control group were also diagnosed with 
both chronic pyelonephritis and chronic cholecystitis. 
The age of the gravidas ranged from 18 to 33 years in 
both groups; the groups were comparable in terms of 
the average age. Statistically significant differences 
in the height and weight of infants between the con- 
trol group and the group of patients were found. The 
groups also differed in terms of arterial blood pressure 
and time of birth. 

Collection of placental samples 

We examined the distal (maternal) portion of the pla- 
centa. The tissues were sampled immediately after de- 
livery (sample ischemia time did not exceed 10 min). 
Placental tissue samples were taken from the central 
areas close to the umbilical cord at a placental depth of 
0.5 cm. The samples were collected from macroscopi- 
cally normal sections of the placenta (without hemor- 
rhage, calcification, necrosis, or fibrin deposition) with- 
out intervening large vessels, washed with saline to 
remove the residual maternal blood and amniotic fluid, 
immediately immersed in RNAlater (Ambion, UK), and 
transferred to be stored at ~80°C until the RNA iso- 
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Table 1 . Characteristics of the examined groups 



Parameters 


PE, N = 10 


Control group, 
iV= 11 


p- value* 


Mean age, years 


ZD = Z 


Zo — o 


n o/i i 


Mean weight, kg 


en _i_ n 

OU ± I 


co -i- a 
OZ = 0 


U.oz4 


Body mass index, BMI 


Zo ± 4 


OQ -1- Q 
ZO 3= O 


n QQO 

U.ooZ 


The mean maximum systolic blood pressure, mm Hg 


1 £9 -1-10 

ioz iy 


101 -1- *3 

1Z1 — 6 


n nnn i 
U.UUU1 


The mean maximum diastolic blood pressure, mm Hg 


lUrt — 1 O 


fin -i- 4 
ou — 4 


U.UUU 1 


Delivery time, weeks 


38 ± 1 


40 ± 2 


0.009 


Birth weight, g 


2783 ± 560 


3549 ± 345 


0.004 


Length at birth, cm 


50 ± 4 


53 ±2 


0.021 


Premature birth, % 


50 


0 


0.012 


Chronic diseases, % 


60 


50 


0.575 



* The significance level was determined by comparing the groups using the Mann-Whitney test or Fisher's exact test. 



lation procedure. A histological examination revealed 
chorionic villi and decidua tissue with fibrinoid necro- 
sis foci and small calcifications in all biopsy specimens 
(Fig. 1). 

RNA Isolation 

Tissue samples (100-200 mg) were homogenized using 
TissueLyser (Qiagen) in Trizol; RNA was then isolated 
using the standard protocol. The concentration of total 
RNA was determined using Nanodrop ND-1000 based 
on absorbance at 260 nm in water. The quality of the 
samples was monitored using an Agilent 2100 Bioana- 
lyzer capillary electrophoresis system (Agilent Tech- 
nologies Inc., Palo Alto, USA) and spectrophotometric 
scanning. 

Microarray analysis 

A genome-wide profile of gene expression in the pla- 
cental tissue was determined using hybridization on 
HT-12 BeadChip microarrays (Illumina, USA) contain- 
ing information about more than 48,000 transcripts. Af- 
ter hybridization, the microarrays were scanned on an 
Illumina BeadArray Reader device. The raw data were 
converted into mean values of the signal intensity for 
each sample (Sample Probe Profile) using the BeadS- 
tudio v3 software package (Illumina). 

Bioinf ormatic analysis 

The data were analyzed in an R software environment 
using the limma program package [12]. Nonparametric 
background correction followed by quantile normaliza- 
tion (neqc function) was performed for the entire data 
set. The specimens that were identified in all the sam- 




Fig. 1 . Micrograph of one of the studied placental biopsy 
samples. Hematoxylin- and eosin-stained 



pies of at least one of the experimental groups (detection 
p-value <0.01) were further considered. A differential 
expression analysis was performed using multiple linear 
regression and moderated t-statistics [12], including the 
assessment of the quality weights of microarray reading 
[13] and Benjamin-Hochberg multiple testing correction 
(FDR) procedure. A 1.5-fold (or greater) change in the 
level of gene expression (FC - fold change) was consid- 
ered to be significant at the adjusted significance level 
of p < 0.1. Functional annotation and functional cluster 
analysis of the groups of differentially expressed genes 
(DEGs) were performed using the DAVID (Database 
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Fig. 2. Heatmap of DEGs (FDR < 0. 1 ; FC < 1 .5). Each column represents a sample; each row represents DEG. Samples 
from PE patients are shown in pink; samples from the females of the control group are shown in blue. The color scale of 
the heatmap indicates the deviation of the normalized expression level in the cell from the mean value for the row 



for Annotation, Visualization and Integrated Discovery) 
web-based tool with the standard values of clustering 
parameters and enrichment score EASE s 0.01 [14]. Con- 
struction of gene networks was performed using the 
STRING 9.0 program (Search Tool for the Retrieval of 
Interacting Genes) [15]. 

This study was approved by the Ethics Committee 
at the Research Institute of Medical Genetics, Siberian 
Branch of the Russian Academy of Medical Sciences. 



RESULTS 

Our analysis identified 63 genes with significantly dif- 
ferent expressions (FDR <0.1; FC a 1.5) in the placen- 
tal tissue of females with PE and physiological preg- 
nancy (50 DEGs with an increased expression level and 
13 DEGs with a decreased expression level). The DEG 
cluster, whose expression was increased in PE, includes 
not only known candidate genes that have previously 
been identified in many genome-wide studies of the ex- 
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pression profiles of the placental genes in preeclamp- 
sia (e.g., LEP, BHLHB2, SIGLEC6, RDH13, BCL6), but 
also new potential candidate genes (COR02A, SYDE1, 
PLIN2, CEBPA, HK2, NDRG1, ERRFI1, EFNB1, GFOD2, 
NCOR2, HMHA1, HERPUD1, KIF2A), whose association 
with the development of PE has been established either 
in a few studies [16-21], or was done so for the first time 
in our work. The products of some of these genes, based 
on current knowledge on their functional features, can 
be involved in the etiopathogenesis of PE. 

Figure 2 shows the heat map with the results of a 
hierarchical clustering of females according to the ex- 
pression level of 63 DEGs. It can be seen that all PE pa- 
tients but one fall into one cluster, while females with 
physiological pregnancy fall into the other one. One PE 
sample was assigned to the control group probably due 
to the significant interindividual variability of the tran- 
scription levels of the placental tissue genes. A similar 
phenomenon was observed in several human cell lines: 
in particular, in the cell lines of the hepatocyte tran- 
scriptome [22, 23]. 

Table 2 shows the data related to the most significant 
DEGs (FC > 2, FDR <0.01). The presence of several genes 
whose products are involved in the transcriptional regu- 
lation (BHLHB2, ZNF175, ANKRD37, BCL6) in this list, 
as well as a significant increase in the expression level of 
the LEP gene and its receptor gene SIGLEC6 during the 
development of PE, is of special interest. 



We analyzed DEG using the DAVID online resource 
to study the biological processes associated with the 
development of PE (Fig. 3). The major categories of 
molecular functions of the protein products of these 
genes include responding to various stimuli, immune 
processes, regulation of cell communication, intracel- 
lular signaling cascades, etc. The analysis of metabolic 
pathways including DEGs has shown that cytotoxicity 
pathways associated with NK-cells, transendothelial 
migration of leukocytes, and signaling pathways medi- 
ated by GTPase activators are probably involved in the 
molecular mechanisms of PE. 

Protein-protein interactions of DEG products were 
analyzed to identify the possible relationships with 
DEGs (Fig. 4). The associations in the constructed net- 
work are mainly based on "text mining" (mentioned 
in the abstract of one article). The coexpression clus- 
ter that includes the RAC2, CYBA, TYROBP, HMHA1, 
ITGB2, LYN and LCP1 genes should be mentioned. In 
addition, LEP and its receptor SIGLEC6 and ephrin 
with its kinase LYN are of certain interest. 

Our study also revealed features of differential gene 
expression depending on PE severity (Table 3). A to- 
tal of eight DEGs were identified (FDR <0.1; FC =s 1.5), 
whose expression levels were significantly different in 
moderate and severe forms of the disease. In our opin- 
ion, HSPA1A encoding the highly conserved heat shock 
protein 70 (HSP70) and BAG3 encoding Bis, a Bcl-2 
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Fig. 4. Protein-protein 
interactions between 
DEG products. The 
proteins are shown as 
circles; the color line 
between these circles 
indicates the evidential 
category of protein-pro- 
tein interaction: yellow — 
literature data ("text min- 
ing"), black -according 
to the analysis of gene 
coexpression, purple — 
the experimental results, 
blue - evidence from the 
databases, pink — cumu- 
lative evidence 




Table 3 . List of differentially expressed genes (FDR <0. 1 ; FC =s 1 . 5) in moderate to severe preeclampsia 



Gene 


Change in gene 
expression level 


FC 


FDR 


Gene product 


HSPA1A 


t 


6.44 


0.079549 


Heat shock protein 70, HSP70-1A 


BAG3 


T 


2.14 


0.073131 


Bcl-2-associated athanogene 3 


SNHG8 


t 


1.78 


0.04105 


Small nucleolar RNA 8 


LOC729660 


I 


2.63 


0.010437 


No data 


LOC728457 


i 


2.43 


0.010437 


No data 


AP0C1 


I 


2.28 


0.04433 


Apolipoprotein CI 


LOC401357 


i 


2.27 


0.010399 


No data 


LOC100128326 


I 


1.92 


0.079549 


No data 



binding protein, are of the greatest interest. The main 
function of the Bis protein is inhibiting the chaperone 
activity of the HSP70/HSC70 complex. 

A comparative analysis of gene expression profiles in 
the placental tissue of females with moderate PE and in 
the control group revealed 56 transcripts of 52 genes, 
whose transcription levels differ significantly in these 
populations. Changes in the expression profile were 
more pronounced in severe PE: a significant increase in 
the expression of 55 genes and a decrease in the expres- 
sion of 35 genes compared to physiological pregnancy 
were observed (Fig. 5). It should be noted that, along 
with a small amount of common genes (21 genes) that 
were differentially expressed in both severe and mod- 
erate PE, more than 60 DEGs were specific only to the 



severe form of the pathology. The results of a functional 
annotation of these genes in the DAVID web-resource 
point to a number of biological processes that are sta- 
tistically significantly associated with the development 
of severe PE, such as processing and presentation of 
peptide or polysaccharide antigens and protein folding 
(Table 4). An analysis of the metabolic pathways that 
involve these genes also demonstrates the important 
role of the mechanisms of processing and presentation 
of antigens in the molecular pathogenesis of severe PE 
(according to the KEGG and BIOCARTA databases). 

DISCUSSION 

The placenta is the key in understanding the physi- 
ological processes associated with pregnancy. It is im- 
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Table 4. Main biological processes that involve differentially expressed genes characteristic of severe preeclampsia 



Categories of biological processes 


Gene 


p* 


Processing and presentation of peptides or polysaccharide 
antigens via class li iVLtiC molecules (IjU.00U<sd04j 


HLA-DPA1, CD74, HLA-DMA, HLA-DRA 


0.0421 


Processing and presentation of exogenous peptide 
antigens (GO:0002478) 


HLA-DMA, CD74, HLA-DRA 


0.0453 


Chaperone-mediated protein folding (GO:0051085) 


ER01L, HLA-DMA, CD74 


0.0467 


De novo post-translational protein folding (GO:0051084) 


EROIL, HLA-DMA, CD74 


0.0478 


Reactions of unfolded protein molecules (GO:0006986) 


EROIL, HSPH1.HSPA1A, HERPUD1 


0.0489 



* Significance level including Benjamin-Hochberg multiple testing correction, which characterized the accuracy of the as- 
signment of this set of genes to a certain biological process. 



portant to characterize the genes essential for placen- 
tal function to understand the mechanisms underlying 
normal and pathological gestation. The results of this 
work show that the identified DEGs belong to several 
biological processes associated with immune responses, 
cell-cell interactions, and responses to various stimuli. 
It should be noted that the analysis performed using 
the module for functional annotation clustering of the 
DAVID bioinformatical resource made it possible to 
identify 16 clusters. However, only one of them had an 
enrichment score of over 2. This cluster includes five 
genes {KRT19, RAC2, LIMCH1, BCL6, LCP1) involved 
in the biological processes related to the organization 
of the actin cytoskeleton (GO: 0030036; GO: 0030029; 
GO: 0007010). Studying the functional role of the ac- 
tin cytoskeleton is one of the important trends in the 
study of cellular signaling mechanisms. Numerous ex- 
perimental data published over the past few years pro- 
vide evidence of the fact that actin is involved in the 
regulation of gene expression and mediates it by par- 
ticipating in transcription elongation, assembly of the 
preinitiation complex, mRNA maturation and export, 
chromatin reorganization, and other processes [28, 29]. 
In this context, the increase in the expression level of 
the COR02A gene seems interesting. The product of 
this gene, coronin 2A, belongs to the family of actin- 
binding proteins and mediates the Coro2A/actin-de- 
pendent mechanism of derepression of the inflamma- 
tory response genes [27]. 

We found no association between the development of 
PE and such canonical pathways as abnormal apoptosis 
and angiogenesis as described in several papers [16, 19, 
30, 31]. This is probably attributable to the interethnic 
variability of the gene expression profiles in the placen- 
tal tissue due to the population differentiation of the 
regulatory regions of the genome, or due to the differ- 



Severe preeclampsia 




Fig. 5. Venn diagram showing the results of the gene 
expression profiling in moderate and severe preeclamp- 
sia and in physiological pregnancy. DEGs - genes that 
are differently expressed in females with preeclampsia 
and physiological pregnancy (control group). The arrow 
shows the increase (T) or decrease (A) in gene expression 



ent criteria (population size, delivery time, severity of 
the disease, etc.) used in the formation of the examined 
groups. The different placental localization of the biop- 
sy samples used in individual studies of the transcrip- 
tome in PE is another factor that apparently affects 
the occurrence of these contradictions. Thus, high- 
throughput sequencing (RNA-Seq) revealed signifi- 
cant differences in the gene expression profiles in the 
amnion, chorion, and decidua of the human placenta 
[32]. Similar findings were previously arrived at when 
performing a microarray analysis of the transcriptome 
patterns in different portions of the placenta [33]. 

Despite the aforementioned differences in the re- 
sults of the DEG functional annotation, it remains of 
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interest that changes in the expression levels of some 
DEGs identified in our work were also described in 
other studies (Table 5). 

Thus, a significant increase in LEP gene expression 
in preeclampsia was observed in almost all genome- 
wide studies of gene expression profiles in the placen- 
tal tissue. Leptin, the product of this gene, is one of 
the new serum markers of PE. It is known that leptin 
belongs to adipocyte-specific cytokines that regulate 
energy metabolism and are involved in various meta- 
bolic and neuroendocrine processes [39]. The studied 
group of PE patients did not differ significantly from 
the control group of patients in terms of weight and 
body mass index. None of the patients had abnormal 
weight gain during pregnancy, and therefore it can 
be assumed that the contribution of leptin to the de- 
velopment of PE is determined by other functions of 
this protein. Placental leptin is known to be involved 
in providing the flow of nutrients to the fetoplacental 
complex and to induce the trophoblast proliferation by 
inhibiting apoptosis [40, 41]. Thus, an increase in the 
leptin level in the placenta may be a compensatory 
mechanism directed against endothelial dysfunction, 
which is observed in PE. Meanwhile, it was shown that 
leptin is involved in the activation of the sympathoad- 
renal system, which contributes to arterial hyper- 
tension, the main symptom of PE [42]. In addition, an 
important immunomodulatory function of leptin was 
found, which may also contribute to the development 
of pregnancy failure [43]. Despite the intensive studies 
devoted to LEP gene expression, only some works have 
focused on the analysis of the hereditary variability of 
this gene and its role in changing the transcription level 
and the structure of susceptibility to pregnancy failure. 
It was shown that carriers of the AA genotype of the 
rs2167270 locus (G19A) located in the promoter region 
of the LEP gene have elevated levels of expression of 
this gene in the blood, as well as the risk of PE and hy- 
pertension [44]. Association between another polymor- 
phism (G2548A) localized at the LEP gene promoter 
with gestational diabetes [45] was found in the Czech 
population. Along with this, several studies [37, 46, 47] 
have revealed significant hypomethylation of this lo- 
cus, as well as dysregulation of the placental epigenome 
during the development of PE. 

The increase in the expression of the long non-coding 
RNA 284 gene in the placental tissue of PE patients, 
which was observed in our work, is of particular inter- 
est in the context of the role of epigenetic dysregula- 
tion in the formation of this pathology. It has recently 
been shown that long ncRNAs perform vital regulatory 
functions in cells. In particular, it is assumed that they 
can function as a module scaffold in the specific, highly 
ordered organization of ribonucleoprotein complexes 



and induction of epigenetic changes in these loci. Some 
long ncRNAs can bind chromatin using remodeling 
enzymes and then participate in the local chromatin 
modification, e.g. in DNA methylation, by initiating or 
repressing transcription. This RNA class can partici- 
pate in the binding of transcription factors and inhibit- 
ing gene expression [48]. 

The association between overexpression of the 
BAG3 and HSPA1A genes and the development of se- 
vere PE seems to be of interest. The protein product 
of the BAG3 gene is known to compete with Hip co- 
chaperone for binding to the ATPase domain of the 
HSP70/HSC70 complex and thus inhibits the chaper- 
one activity of the heat shock protein 70 (Hsp70), the 
product of the HSPA1A gene, whose expression is sig- 
nificantly elevated in severe PE (more than sixfold as 
compared to moderate PE and eightfold as compared 
to the control population). It is known that Hsp70 per- 
forms various functions. Improving the resistance of 
the protein biosynthesis apparatus to damaging in- 
fluences and chaperone activity are the most signifi- 
cant of them. In addition, there is data indicating that 
Hsp70 participates in protein transport, conduction 
of the intracellular signal, and protease-dependent 
degradation [49]. It should be noted that according 
to the results of the functional annotation of DEGs, 
processing and presentation of peptide and polysac- 
charide antigens and chaperone-mediated protein 
folding are the principal processes characteristic of 
severe PE. Since the heat shock protein Hsp70 is ca- 
pable of forming complexes with non-folded proteins 
and a wide variety of peptide fragments that are pre- 
cursors of the antigenic peptides presented on the cell 
membrane along with other class I and II MHC mol- 
ecules [50], it is reasonable to assume that the immu- 
nological control mechanisms of trophoblast invasion 
in the uterine wall and the immunological tolerance 
factors in the mother-fetus system play a key role in 
the pathogenesis of severe PE. The pathological ef- 
fect of these factors can lead to gestational complica- 
tions. Furthermore, many heat shock proteins exhibit 
immunoregulatory activity, stimulate the maturation 
of dendritic cells, and induce some proinflammatory 
cytokines [51]. These properties of these proteins may 
also contribute to the mechanisms of severe PE. 

The statistically significant decrease in the expres- 
sion of the APOC1 gene in severe PE, which was re- 
vealed in the present study, is probably due to the de- 
velopment of oxidative stress in the blood vessels of the 
placenta or the recently discovered immunosuppres- 
sive properties of the CI apolipoprotein encoded by this 
gene [52]. It was previously shown that the serum of PE 
patients has a high level of triglyceride-rich lipopro- 
teins, which can promote endothelial dysfunction [53, 
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Table 5. Differentially expressed genes identified in this study whose association with preeclampsia has been previously 
shown in studies focused on the placental tissue transcriptome 





Gene 


Gene product 


FC 


Significance 
level 


T^.tVinin nnni ila tinTit; 


TJpfpypnr'p 


1 


LEP 


Leptin 


10.94 


< 0.0001* 


Japanese 


[30] 


8.58 


0.036* 


Chinese 


[34] 


40.11 


1.35 x 10" 9 


Caucasian 


[16] 


5.52 


0.0020 


Caucasian 
Afro- American 
JMongoloid 
Spaniard 


[18] 


108.9 


< 0.0001 


Caucasian 


[35] 


4.4 


< 0.0001 


Korean 


[36] 


> 1.5 


< 0.05 


Chinese 


[37] 


11.79 


<0.01* 


American 


[38] 


2 


BCL6 


Transcription factor - zinc finger 
protein 51 


1.78 


0.0154 


Caucasian 

Afro- American 

Mongoloid 

■ i 

Spaniard 


rial 


2.02 


0.0024 


Japanese 


roni 

[30] 


2.24 


O R O / 1 fl-5 

3.00 x 10 


Caucasian 


n 

[16] 


2.60 


< 0.01* 


American 


[38] 


> 1.5 


< 0.05 


Chinese 


[37] 


3 


SIGLEC6 


Immunoglobulin-like lectin 6, which 
binds to sialic acid 




0.02* 


American 


roil 

[31] 


O 1 o 


U.U01 


Caucasian 


n ci 
[16] 


Z.IO 


x n [11* 

< U.01 


American 


rooi 

[38] 


> 1.0 


< u.uo 


Chinese 




A E 

4.5 


u.oiy 


Caucasian 


[35] 


A 
<± 


p n rj j o 


Retinol dehydrogenase 13 




.ISO X 1(J 


American 


TQ 1 1 

[31J 


i.y i 


<- n ni* 
^ U.U1 


American 


rooi 
[doj 


> 1.5 


< 0.05 


Chinese 


[37] 


5 


NDRGl 


Cytoplasmic protein belonging to the 
hydrolase superf amily 


2.02 


0.0001 


Japanese 


[30] 


2.67 


1.12 x io- 5 


Caucasian 


[16] 


6 


BHLHE40 


Transcription factor with helix-loop- 
helix domain 


1.95 


0.0004 


Japanese 


[30] 


3.08 


2.18 x 10" 5 


Caucasian 


[16] 


7 


KRT19 


Keratin 19 


1.75 


0.0071 


Japanese 


[30] 


2.28 


1.59 x IO" 5 


Caucasian 


[16] 


8 


GPT2 


Alanine aminotransferase 2 


2.45 


3.70 x 10" 5 


Caucasian 


[16] 


9 


PPP1R12C 


12 A regulatory subunit of phos- 
phatase 1 




2.16 x 10 8 * 


American 


[31] 


10 


CEBPA 


CCAAT/enhancer-binding protein a 




2.52 x 10- 8 * 


American 


[31] 


11 


HK2 


Type 2 hexokinase 


3.90 


3.87 x 10" 6 


Caucasian 


[16] 


12 


HMHA1 


Minor histocompatibility antigen HA1 




1.23 x IO" 8 * 


American 


[31] 


13 


PVRL4 


Nectin 4 


2.54 


3.62 x io- 5 


Caucasian 


[16] 


14 


SASH1 


SAM- and SH3-domain-containing 
protein 1 


2.54 


1.22 x 10- 7 


Caucasian 


[16] 


15 


SH3PXD2A 


SH3- and PX-domain-containing 
protein 2A 


>1.5 


< 0.05 


Chinese 


[37] 


16 


SYDE1 


GTPase activating protein, homolog 1 


1.55 


< 0.01* 


American 


[38] 



* Significance level including multiple testing correction. 
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54]. Meanwhile, the blood level of E and Al apolipopro- 
teins is decreased in this pathology [55, 56]. In addition, 
a protective role of the APOE e2 allele in the Kurd pop- 
ulation was demonstrated, which is related to the high 
antioxidant properties of this allele according to the au- 
thors [57]. We failed to find any information about an 
association between APOC1 polymorphisms and ges- 
tational complications. However, it was shown that the 
insertion-deletion polymorphism at —317 position of the 
promoter region of this gene (rsl 1568822) is associated 
with Alzheimer's disease, while the rs4803770 marker 
is associated with the coronary heart disease [58, 59]. 
Since the APOC1 gene localizes in the same cluster as 
the APOE gene, it is assumed that these associations 
are due to the strong linkage disequilibrium between 
these genes [60]. However, we found no statistically sig- 
nificant changes in the expression level of the APOE 
gene. Therefore, it is reasonable to consider the APOC1 
gene to be an "independent" new candidate gene for 
PE. However, this assumption needs confirmation. 

Thus, the findings indirectly confirm the immunolog- 
ical hypothesis of the development of severe PE, which 
postulates the key role of immune competent cells (B 
lymphocytes, monocytes, dendritic and NK cells) in the 
pathophysiology of this disease. This theory assumes 
that the etiopathogenesis of PE is triggered by insuffi- 
cient trophoblast invasion into maternal spiral arteries, 
which is associated with a decreased expression of HLA 
antigens and "aggression" against NK cells. This results 
in reduced placental perfusion and development of hy- 
poxia at the mother/fetus boundary, which, in turn, 
triggers the activity of pro-inflammatory cytokines, 
leading to endothelial dysfunction [61]. B cells may also 
contribute to the development of preeclampsia by pro- 
ducing anti-adrenoceptor antibodies. 

CONCLUSIONS 

The present work is the first Russian genome-wide 
study of differential gene expression in the placental tis- 
sues in normal and complicated pregnancies. The results 
indicate that some processes can play an important role 
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